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Abstract 

The yrast states of transitional nuclei are described as quadrupole waves run- 
ning over the nuclear surface, which we call tidal waves. In contrast to a rotor, 
which generates angular momentum by increasing the angular velocity at ap- 
proximately constant deformation, a tidal wave generates angular momentum 
by increasing the deformation at approximately constant angular velocity. The 
properties of the tidal waves are calculated by means of the cranking model 
in a microscopic way. The calculated energies and E2 transition probabili- 
ties of the yrast states in the transitional nuclides with Z— 44, 46, 48 and 
N = 56, 58, 66 reproduce the experiment in detail. The nonlinear response 
of the nucleonic orbitals results in a strong coupling between shape and single 
particle degrees of freedom. 



Most of the theoretical interpretation of the low-spin structure of nuclei in 
the transitional region between spherical and well deformed shape has been 
carried out in the framework of phenomenological models based on the Bohr - 
Hamiltonian[l|,[3,[3| or using the algebraic IBM concept 0, M, [l] ■ Although these 
studies account for many low-spin aspects of the collective quadrupole degree of 
freedom in an impressive way, the connection with the underlying microscopic 
structure has remained a challenge. On the other hand, the various versions 
of mean field theory are very successful in predicting shell closures, static de- 
formations, and delineating the regions of spherical and deformed shapes (for 
recent reviews cf. e.g. @, UQ)- Yet the large amplitude collective dynamics of 
transitional nuclei starting from the microscopic mean field remains a challenge 
to nuclear theory. At present there are two mature approaches of this type: the 
combination of mean field potential energy with cranking mass parameters (see 
review [l(J) and the Generator Coordinate Method (see U, 12] for recent devel- 



opment). Both approaches use the adiabatic approximation that the collective 
quadrupole excitations do not mix with the quasiparticle excitations, which is 
often not justified. 

In this Letter, we demonstrate that the yrast states of transitional nuclei 
have the very simple structure of a running surface wave (tidal wave), which 
permits us to calculate their properties microscopically by means of the cranked 
mean field theory without resorting to the adiabatic approximation. First we 
introduce the concept of tidal waves for the collective quadrupole vibrations of 
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Figure 1: Potentials of a vibrator (V), soft rotor (R), and transitional nucleus (T). 

a classical droplet with irrotational flow, which is discussed in 0, Q . Then we 
will extend it to the microscopic mean field theory and carry out calculations 
for a group of transitional nuclei around A=110. 

The quadrupole displacement of the surface with respect to the sphere is 
R(t) — R sp h = Ylu a n(t)Y2)i(^, The displacement is expressed by the familiar 
deformation parameters f3, 7 and the vector <fi = n(f> (rotation about the axis n 
by the angle (/>), which fixes the orientation of the quadrupole shape. 



^ = u% cos 7 + (yiS) + ul_S)) 
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functions.) The kinetic energy is 
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J i =4^/3 2 sin 2 ( 7 -i^) 
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and the potential energy V{(3, 7) depends only on the deformation parameters. 
The equations of motion are derived from the Lagrangian T — V in the standard 



way. The solution with maximal angular momentum J in z-dircction for given 
energy (yrast mode) is a wave that runs over the spherical surface with the 
constant angular velocity uj. In the frame of reference rotating with the velocity 
ui the surface has a static deformation. The deformation parameters /?, 7 are 
the equilibrium values at the minimum of the energy 

EW,l) = + (4) 

J = 4Df3 2 sin 2 (~f~ 2-), (5) 

where the angular momentum J = Juj and uj — dE/dJ. 

The location of the droplet surface in spherical coordinates is given by 

R(t) = R a [l + V2/3shxjcoa(24>-2u}t)Y 2 2('&,<p = 0). (6) 

The inset of Fig. [2] illustrates this mode, which looks like the tidal waves on 
the oceans of our planet caused by the gravitational pull of moon and sun. In 
allusion, we suggest using this name for the running waves. Concerning the 
motion of the surface, it does not differ from rigid rotation. The difference 
between a rotor and a tidal wave is the way angular momentum is generated. 
The rigid rotor has a fixed deformation and a constant moment of inertia J . 
The angular momentum J — Juj is generated by increasing the angular velocity 
uj. The tidal wave runs with a fixed angular velocity u), while the angular 
momentum is gained by increasing the moment of inertia by changing /3, 7. 
Of course, real nuclei lay between these idealized cases. Thus, we suggest the 
following classification: A rotor generates most of its angular momentum J by 
increasing the angular frequency uj, while the moment of inertia J does not 
change much. A tidal wave generates most of J by increasing J ', while uj does 
not change much. 

Let us discuss some cases in more detail. 
-Harmonic vibrator (V in Figs. I1I2|) : 

V = ^V\ (7) 

Minimizing the energy one finds 

T, /3 2 = -4=, J = 4Df3 2 = ^=, (S) 



' 6 

uj = - = -J-, E = ujJ = n- = Cf3 2 . (!» 

The wave travels with an angular velocity being one half of the oscillator fre- 
quency ft. The angular momentum is generated by only increasing the moment 
of inertia, which is a linear function of J. This is the ideal tidal wave. The 
solution with the same energy but J = is R(t) = Ro(l + cos(r2i)Y~2o($)), 
which is the familiar oscillation between prolate and oblate shapes. 
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-Rigid axial rotor (R in Fig. [5]): 




V 9 <5(cos(3 7 ) - 1) (10) 

The potential has a minimum at j3o, where V ~ C((3 — (3o) 2 /2. Assuming C 
being infinite large one finds 

7 = -2tt/3, P = A), Jo - 3D/3j| = const. (11) 

-The soft rotor: (R in Fig. [TJ 

With a finite C the minimization gives 

2 7 2 

/3=^o(l + T? ^ + 0(J 4 )), (12) 

o r2 

J = Ml+ JoCW +0{j4)) ' (13) 

The moment of inertia is a slowly increasing quadratic function of J. 
-Transitional nucleus (T in Figs. Il|2|) : 

The yrast energies are well approximated by the expression suggested [3] , which 
gives (J= y/WTTj) 

F(n J 2 7 (Qo + exJ) 2 

"V' (2(60 + 6x7)'^ (60 + 6! J/2) l±4J 

As expected for a transitional nucleus, the angular momentum is gained by 
increasing both J and w. From a vibrational perspective, the increase of ui 
reflects the anharmonicities of the motion, from the rotational perspective, the 
increase of J reflects the softness of the rotor. Assuming the energies are derived 
from minimizing (j4|) and axial shape (7 = 0) one finds 

P(J) = 6 + 61 J 

^5(60 + 61 J/2)' 
T/ . n J 2 (60 + 6! J/2) 

F(J) = (60 + 6, J) 2 ' (16) 

which is a parametric form for V(/3/(3o). It is displayed in Fig. [T]as T, which 
has a form intermediate the vibrator (V) and rotor (R). Using Qt — Q/3//3q 
for the transition quadrupole moment of E2-transitions, one may calculate its 
spin dependence by means of (fT5)l . Tab. [T] shows the case of 154 Gd, which has 
been classified as a transitional point nucleus [l6j]. The parameters 60 = 
21.5MeV~ 1 and 6i = 1.18MeV _1 place the nuclide slightly above the T -line 
in Fig . [3J The values for Q t are very similar to the ones of the X(5) limit given 
in [16|, which is consistent with the the shape of the potential V(/3//3o)- The 
increase of Qt is larger than that of the experimental values Qt,exp> which points 
to degrees of freedom other than the quadrupole deformation contributing to 
the growth of J{I)- 
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Table 1: Energies E (in MeV) and transition quadrupole moments Q (in eb) for 154 Gd. The 
subscript pm is used for the phenomenological model given by j 1 4[ ) and II15I I and e for the 
experiment 15, 16]. The columns /3//3q and V display the deformation potential in parametric 
form, 
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The transition between rotation and tidal waves is gradual. Let us classify 
the yrast mode of nuclei as an (anharmonic) tidal wave when most of the angular 
momentum is generated by the increase of the moment of inertia and as (soft) 
rotation when most of the angular momentum is generated by the increase of the 
angular velocity. Experimental functions J(w) and J{ J) can be derived from the 
level energies by the familiar procedure indicated in the caption of Fig. [2j which 
shows the Nd-isotopes. For the well deformed N = 94 isotope, the increase of 
J remains moderate compared with J7(0), and it grows approximately with 
J 2 , which is characteristic for a rotor. The frequency ui starts with a small 
value and grows steadily. The isotopes with N = 84 — 88 belong to the region 
of tidal waves. The frequency w starts with a large value and stops growing, 
while J starts with a small value and increases substantially. As expected, 
N = 84 comes closest to the spherical vibrator limit. The experimental functions 
J(uj) and J(J) become increasingly irregular when approaching the shell closure 
at N — 82, which indicates that the the motion looses collectivity becoming 
progressively entangled with the single particle degrees of freedom. The limit of 
a good vibrator is not realized in nuclei. The isotop N — 90 marks the border 
between a tidal wave and a soft rotor, which is in accordance with the onset of 
a substantial ground state deformation. 

The discussion of the liquid drop model for the quadrupole surface modes 
suggests that one may calculate the yrast states in a microscopic way by finding 
the selfconsistent mean field solution in a rotating frame of reference, i. e. one 
may use the cranking model for calculating the tidal wave modes in transitional 
and vibrational nuclei. In fact, pjj demonstrated that the RPA equations for 
harmonic quadrupole vibrations of nuclei with a spherical ground state are a 
solution to the selfconsistent cranking model if the deformation of the mean 
field is treated as a small perturbation. In this paper we solve numerically the 
cranking problem without the small deformation approximation, which allows 
us to describe the yrast states of all nuclei in the range between the harmonic 
vibrators and rigid rotors. Such microscopic approach takes into account not 
only the quadrupole shape degrees of freedom, as in the above discussed phe- 
nomenological model, but also the the quasiparticle degrees of freedom, which 
play an important role to be discussed below. Although we only study even-even 
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nuclei in this Letter, we want to point out that the method can be applied to 
odd-A and odd-odd nuclei without any problems. 

We use the Shell Correction version of the Tilted Axis Cranking model (SC- 
TAC) in order to study the tidal waves in the nuclides with Z =44, 46, 48 and 
N =56, 58, 60, 62, 64, 66. The details and parameters of SCTAC are described 
in [l8| . For the yrast states of even-even nuclei, the axis of rotation coincides 
with a principal axis (x). The condition < J x >= J = I is used to fix u, which 
is solved numerically for each point on a grid in the £2 — 7 plane. In solving 



it is crucial to employ the diabatic tracing technique described in [18| , which is 
illustrated in Fig. SJ The ground state (g-) configuration indicated by the full 
dots is followed through crossings of the hn/2 routhians. The avoided crossings 
between the positve parity routhians around Tuuj — 0.53 MeV are taken into 
account adiabatically: The occupation is kept on the lower of the two levels. 
The yrast states with high spin correspond to the s-configuration indicated by 
the open circles 

The interpolated function _E(e2,7, J) is minimized with respect to deforma- 
tion parameters for fixed angular momentum. The final value of w(J) for the 
equilibrium deformation is found by interpolation between the values at the grid 
points. It is shown as the theoretical points in the Figure [5] It is noted, that 
the standard technique of finding the selfconsistent solutions by iteration at 
fixed uj becomes problematic for tidal waves close to the harmonic limit because 
the total Routhian is nearly deformation independent. The pairing parame- 
ters are kept fixed to A„ = A p = 1.2 MeV, and the chemical potentials are 
A adjusted to the particle numbers at ui = 0. The B(E2, 1— > I — 2) values 
are calculated for J = 7 as described in [18] for 9 = 90° and multiplied by 
7(7 — l)/(7 — 3/2)(J— 1/2). The factor corrects for the low-spin deviation of the 
B(E2) values from the high-spin limit used in the TAC code. It represents the 
Clebsh-Gordan coefficient appearing for axial K — bands [l|, which is justified 
by the fact that for most nuclei 7 < 10°. We slightly correct 

J{uj) = Jsctac(u) + WOMeV-^l sin 2 ( 7 - tt/3), (17) 

which adds less than 0.5 Ti to the angular momentum. About half of it takes into 
account the coupling between the oscillator shells and another half is expected 
to come from quadrupole pairing, both being neglected in SCTAC. The changes 
due to corrections are smaller than the radius of the open circles in Figs. [5] and 

El 

The ground state energies are shown in Fig. [3] They are calculated using 
the modified oscillator potential with the parameters given in Ref. Q3 • The 
tendency to deformation increases with the numbers of valence proton holes 
and valence neutron particles in the Z = N = 50 shell. The Z — 48 isotopes 
are spherical, becoming softer with increasing N. The Z = 44 isotopes are 
slightly deformed, very soft, becoming less soft with N. The Z — 46 isotopes 
are intermediate. Figs. [5] and [6] compare the calculated values of J{uj) and 
B(E2, 7— > I — 2) with experiment. The detailed agreement was achieved 
by adding to ground state deformation energies the small adjustment ke\ with 



6 



Table 2: Rotational frequencies in 110 Cd. The subscript o is used for the calculation without 
the correction of the ground state energies, the subscript c for the one with the correction, and 
the subscript e for the experiment. The calculated energy £?(0j)=3.51 MeV, its experimental 
value 1.47 MeV. All energies are in MeV. Data from |15|| 
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k =8.5, 7.9, 8.5, 9.6, 11.9, 14.5 MeV for N= 56, 58, 60, 62, 64, 66, respectively, 
which makes the potentials slightly stiffer. The correction term is substantially 
smaller than the differences between deformation energies calculated from the 
various mean held theories used at present. Tab. [2] compares calculations with 
and without the correction term. 

As seen in Fig. [5J for low-spin (J < 6) the tidal waves in Cd isotopes show 
the linear relation B(E2) oc J that is expected for vibrational nuclei. However, 
the energies deviate strongly from the vibrator limit. For Ru and Pd, the 
B(E2) values start at larger values and increase less with J, such that they do 
not extrapolate to the origin. This reflects the transition to more rotational 
behavior, which becomes stronger with increasing N. The low-spin properties 
reflect the ground state deformation energies in Fig. [3] 

At high spin, there is an irregularity in the J(w) curves, which is caused 
by the rotational alignment of a pair of hu/2 neutrons (change from the g- 
conflguration to the s- configuration in Fig. |4j. It is well known that the 
sharpness of the alignment depends sensitively on the position of the neutron 
chemical potential A„ relative to the hu/2 levels J26|, which causes the famil- 
iar oscillations between rapid alignment (back bending) and more gradual one 
(vertical section) with increasing N. The h±i/2 levels move with the changing 
deformation, and differences of the deformation are the reason why J{u) sharply 
back bends in Z = 48, A = 62 while it gradually grows in Z = 44, N = 62. 
Also normal parity neutron orbitals contribute to the J(uj) in a non-linear way, 
in particular at the avoided crossings between them (around Tiuj = 0.53 MeV 
in Fig. [4j. For example, they are responsible that in 110 Cd the point 1 = 8 
has a lower to than the point 1 = 6. The same is observed in 108 Cd [23J] (not 
shown) . The large frequency encountered at small deformation makes also other 
orbitals than the high-j intruders react in a non linear way to the inertial forces. 
For many nuclides the quasiparticle degrees cannot be treated in a perturbative 
way for / > 6, which means the separation of the collective quadrupole degrees 
of freedom becomes problematic. For the nuclei with the smallest numbers of 
valence particles and holes this entanglement of collective and single particle 
degrees appears already for 1 = 4, which is reflected by the irregular curves 
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J(uj) for Z = 48, N = 56, and Z = 60, N = 84 in Figs. [SJ anddj respectively. 

As seen in Fig. [6j the aligned hn/2 quasi neutrons stabilize the deformation 
at a smaller value than reached before, which is reflected by the decrease of 
the B(E2) values. The deformation increases slower than before the alignment, 
which means that the motion becomes more rotational. In fact, it takes the 
character of Antimagnetic Rotation, discussed in 25j. In some cases, the align- 
ment process is extended over several units of angular momentum, keeping uj 
nearly constant. This means the simple picture of tidal waves, as running at 
nearly constant angular velocity and gaining angular momentum by deforming 
does not account for the full complexity of real nuclei. Other degrees of freedom, 
in particular the single particle ones, may also provide angular momentum such 
that the rotational frequency remains constant. 

102 Pd is a special case. The g- and s-configurations do not mix and are both 
observed up to I — 20. Figs. [5]and|6]show only the g-sequence, which represents 
a tidal wave that gains angular momentum predominantly by increasing the 
deformation, as discussed above for the droplet model. Its amplitude at / = 16 
corresponds to eight stacked phonons! Only for I > 16 the alignment of the 
positive parity neutrons (avoided crossing at Tiuj ~ 0.53 MeV in Fig. [4| becomes 
a comparable source of angular momentum. 

Of course, our cranking approach also allows us to calculate rotational bands 
built on quasiparticle excitations. As an example, Tab. [5] contains the lowest 
two-quasiproton excitation across the Z = 50 shell gap. The excited quasipro- 
tons drive the nucleus to a substantial deformation, which, in contrast to one of 
the yrast tidal wave states, weakly increases with spin. The Ojj" "intruder state" 
has been assigned to this configuration, whose deformed shape coexists with the 
near spherical ground state shape [27j |. The energies of the rotational band are 
well reproduced by the calculation, which however gives a too high energy of 
the 0^" state. 

In conclusion, the yrast states of vibrational and transitional nuclei can 
be understood as tidal waves that run over the nuclear surface, which have a 
static deformation, the wave amplitude, in the frame of reference co-rotating 
with the wave. In contrast to rotors, most of the angular momentum is gained 
by an increase of the deformation while the angular velocity increases only 
weakly. Since tidal waves have a static deformation in the rotating frame, the 
rotating mean field approximation (Cranking model) provides a microscopic 
description. The yrast states of "spherical" or weakly deformed nuclei with 
44 < Z < 48, 56 < N < 66 were calculated up to spin 16. These micro- 
scopic calculations reproduce the energies and electromagnetic E2 transition 
probabilities in detail. The structure of the yrast states turns out to be com- 
plex. The shape degrees of freedom and the single particle degrees of freedom 
are intimately interwoven. The non-linear response to the inertial forces of the 
individual orbitals at the Fermi surface determines the way the angular mo- 
mentum is generated. This lack of separation of the collective motion from the 
single particle becomes increasingly important with decreasing number of va- 
lence nucleons. It can set as low as / = 6 already. Purely collective models 
like IBA |4| or the Dynamical Collective Model [2[ miss this explicit coupling 
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between single particle and collective freedoms. 
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Figure 2: Angular momentum J as function of the angular frequency u> (upper panel) and 
kinematic moment of inertia J as function of J (lower panel) for the Nd-isotopes. The points 
correspond to = (E(I) ~ E(I - 2))/2 and J(I) = - where E(I) is the 

experimental yrast energy. The dashed lines show the cases discussed in the text. 





Figure 4: The Quasi Neutron Routhians for N = 62 and £2 = 0.13, 7 = 0. Full lines: positive 
parity, dash dotted lines: negative parity {hn^)- The full circles show the g-configuration, 
which corresponds to a tidal wave, the open circles show the s-configuration, which has the 
character of antimagnetic rotation. 
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Figure 5: Angular momentum J as function ofjtj 10 angular frequency ui in the mass 110 region. 
The quantities J and w are obtained as in Fig. [2] from the calculated or experimental yrast 
energies E(I). Full symbols show the experiment, open the calculations. Data from [151 . 



